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Abstract 


This  report  describes  a  new  data  analysis  and  interpretation  algorithm  developed  for  a  precision 
meteorological  sensor  that  has  been  fitted  to  one  of  the  Army  Research  Laboratory  s  (ARL) 
prototype  robotics  vehicles.  The  new  algorithm  provides  low-atmosphere  profile  and  turbulence 
information  from  the  retrieval  and  analysis  of  wind  speeds,  wind  directions,  air  pressure,  and 
temperatures  (1  Hz  data).  The  algorithm  is  at  first  formulated  for  data  collected  over  flat  terrain  at 
some  low  reference  height  above  ground.  We  calculate  the  small  interval  (<  4.5  minute)  averages 
and  standard  deviations  for  the  wind  and  temperature  fields  and  apply  a  simple  log-law  model  to 
extrapolate  wind  velocity  profile  information  upward  through  the  sensor.  Thus  we  initiate  a 
coupling  of  meteorology  to  robotics  for  future  military  applications. 
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1.  Introduction 


The  U.S.  Army  Research  Laboratory  (ARL)  is  developing  several  prototype 
robotics  systems  designed  to  carry  combat  sensors  like  infrared,  video,  and 
acoustic  arrays.  In  addition,  an  all-terrain  robotics  vehicle  is  being  fitted 
with  a  precision  meteorological  sensor  to  retrieve  and  construct  local  wind 
field,  temperature,  and  turbulence  information.  It  is  expected  that  instanta¬ 
neous  meteorological  measurements  (with  the  use  of  atmospheric  computer 
models)  will  provide  important  information  to  the  Army's  soldiers  on  the 
performance  of  future  battlefield  systems  and  communications  (e.g.,  http:/ 
/  www.darpa.mil  /fcs  /  index.html,  03-26-02). 

Calculating  low-atmosphere  profile  and  turbulence  information  from  ro¬ 
botics  systems  logically  extends  this  new  capability  and  provides  a  basis 
for  adding  additional  sensors  on  future  testbed  vehicles.  Therefore,  we  re¬ 
port  here  on  a  small  electronic  meteorological  sensor,  which  will  retrieve 
data  at  some  low  reference  height  above  (at  first)  flat  terrain.  In  addition  we 
develop  an  analysis  algorithm  to  calculate  the  small  interval  (<  4.5  min) 
averages  and  standard  deviations  for  the  wind  and  temperature  fields.  We 
then  apply  a  simple  log-law  model  to  extrapolate  profile  information  (~15 
m)  upward  through  the  sensor.  Improved  or  more  sophisticated  models  for 
the  lower  boundary  layer  can  be  implemented  later  on  for  more  complex 
areas,  all  while  trying  to  determine  the  best  configurations  for  single  and 
multiple  robotic  vehicles  (R.  Meyers,  2001  personal  communication).  At  that 
time,  such  algorithms  will  affect  state-of-the-art  theory  and  analysis,  to  in¬ 
clude  corrections  for  inhomogeneous  and  anisotropic  turbulence  conditions. 
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2.  Electronic  wind  velocity  and  meteorological  sensor 


ARL's  Robotics  Systems  Integration  Team  is  investigating  robotics  designs, 
which  include  the  implementation  of  data  display  and  sensor  systems  as 
part  of  the  Army's  command  and  control  architecture.  Compact  packaging 
is  important.  Therefore,  one  of  the  components  being  evaluated  for  an  all- 
terrain  robotics  systems  platform  is  the  Cossonay^  Model  500-C-M  preci¬ 
sion  thermal  anemometer  and  meteorological  sensor.  The  Model  500-C-M 
sensor  is  compact,  automatic,  and  self-contained  with  no  moving  parts  (fig. 
1).  The  Model  500-C-M  is  a  patented,  solid-state  instrument  that  provides 
continuous  digital  (RS232)  signal  output  and  operates  from  a  12v  DC  power 
supply.  The  meteorological  components,  microprocessor,  and  electronics  are 
enclosed  in  a  sealed  weatherproof  casing.  The  sensor  height,  head  width, 
body  diameter,  and  total  weight  are  255  mm,  90  mm  and  70  mm,  and  1.5  kg, 
respectively.  The  Cossonay  500-C-M  sensor  measures  wind  speeds  from  0 
to  40  m/s  (threshold  0.3  m/s),  wind  direction  over  360°,  air  temperature 
from  -40  to  +60  °C  (resolution  0.1  °C),  and  air  pressure  from  600  to  1200 
mbar. 


Figure  1.  Front  view  of 
the  Cossonay  500-C-M 
digital  thermal 
anemometer  and 
meteorological  sensor. 
The  sensor  shown  here 
has  3  platinum  hot  films 
mounted  on  a  small 
ceramic  tube  installed 
vertically  on  the 
instrument  head  and  is 
surrounded  by  a 
protective  wire  cage. 


^  The  use  of  commercial  or  company  names  with  regard  to  electronic  products  does  not 
constitute  an  endorsement  by  the  U.S.  Army. 
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In  general,  thermal  (or  hot-film)  anemometry  is  based  on  the  microanalysis 
of  convection  currents  as  they  are  affected  by  wind  flow  across  heated  plati¬ 
num  wires  encased  by  very  thin  (millimeter)  glass  films  (e.g.,  Thompson 
and  Hassman,  2001).  The  Model  500-C-M  sensor  shown  here  has  3  plati¬ 
num  hot  films  mounted  on  a  small  ceramic  tube  installed  vertically  on  the 
instrument  head  and  is  surrounded  by  a  protective  wire  cage.  The  hot  film 
elements  are  maintained  at  a  high  constant  temperature  above  ambientl 
The  sensor  operates  then  as  explained  above,  i.e.,  heat  transferred  to  the 
atmosphere  by  convection  will  be  greater  for  the  hot  film  element,  which  is 
at  the  leading  edge  of  the  flow,  than  for  the  element,  which  is  at  the  trailing 
edge.  Therefore,  the  differences  in  temperature  across  the  sensor  elements 
are  affected  by  the  magnitude  and  direction  of  the  wind  and  the  tempera¬ 
ture  difference  between  the  sensor  elements  and  ambient  air  temperature. 
(Note:  Ambient  air  temperature  is  detected  by  a  standard  platinum  resis¬ 
tance  thermometer,  colocated  on  the  instrument  head,  as  shown  in  figure 
1.) 

In  addition,  the  Model  500-C-M  sensor  houses  a  precision,  micro-mechani¬ 
cal  pressure  sensor  that  features  similar  technology  to  that  foxmd  in  auto¬ 
mobiles  which  trigger  protective  airbags  upon  impact.  However,  with  re¬ 
gard  to  vibration  or  shock  to  the  instrument  package  in  mobile  applica¬ 
tions,  it  was  suggested  that  future  sensor  designs  may  include  a  reinforced 
metal  instrument  basket,  like  that  found  on  the  Cossonay  Series  270  sen¬ 
sors.  Cossonay  Series  270  crosswind  sensors  have  been  installed  on  tanks 
and  other  military  ground  vehicles  (R.  Jones,  Personnal,  Communication). 
Finally,  the  Model  500-C-M  sensor  will  maintain  accuracy  with  tilt  up  to  20° 
from  vertical  (Model  500  Operating  Manual,  2001). 


2  Ducharme  et  al.  (1994)  reported  on  a  Thermo  Systems  Inc.^  Teflon  coated  hot-film  sen¬ 
sor  that  operated  at  250°  C  above  ambient  air  temperature.  In  contrast,  Thompson  and 
Hassman  (2001)  report  on  a  newer  Sutron  Instruments  Co.^  twin-hot-film  sensor  that  oper¬ 
ates  at  100°  C  above  ambient  air  temperature.  The  Cossonay  500-C-M  sensor  operates  in  the 
range  100-150°  C. 
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3.  Data  analysis  algorithm 


In  this  section  we  present  a  simple  data  analysis  algorithm  for  calculating 
the  small  interval  (<  4.5  min)  averages  and  standard  deviations  from  mea¬ 
sured  wind  speeds,  wind  directions,  air  pressure,  and  temperatures.  We 
follow  the  example  shown  by  an  earlier  program  code  for  partitioning  the 
raw  data  and  filtering  spurious  noise  (Chintawongvanich,  1993).  This  ele¬ 
ment  of  the  calculation  operates  on  the  data  in  the  following  manner.  We 
designate  the  number  of  points  (e.g.,  <  90)  for  the  "small"  averaging  period, 
not  to  exceed  1.5  minutes  of  data  retrieved  at  1  Hz.  Then  we  determine  the 
width  of  the  "large"  averaging  period,  which  will  consist  of  some  number 
(e.g.,  2  or  3)  of  "small"  averaging  periods.  In  total,  a  "large"  averaging  pe¬ 
riod  will  contain  <  4.5  minutes  of  filtered  data.  Then  we  apply  a  noise-rejec¬ 
tion  (i.e.,  spike  filtering)  algorithm,  which  has  the  following  properties.  Once 

the  "small"  averaging  period  mean  (T)  and  fluctuation  (T')  values  for  each 
meteorological  parameter  are  determined,  then  spurious  noise  is  defined 
as  individual  data  points  having  values  greater  than  or  equal  in  magnitude 


to  3  standard  deviations,  e.g.,  30^,  where 


K-r.-r) 


i=l 


n-1 


(Quality  Con¬ 


trol  Software,  2001).  Here  we  consider  the  parameter  air  temperature  (T-), 


where  f  ^  ^  7)  and  f  =T--T  ,  where  n  is  the  number  of  data  points.  As  a 
(=1 

result,  if  data  spikes  are  detected,  the  values  of  those  data  points  are  set  to 
the  mean  for  that  averaging  period. 

The  next  element  of  the  algorithm  is  the  application  of  a  simple  log-law 
model  to  determine  low-atmosphere  wind  velocity  profile  and  turbulence 
information.  As  described  by  Mimn  (1966),  the  vertical  wind  profile  in  the 
absence  of  buoyancy  (i.e.,  in  adiabatic  conditions)  can  be  expressed  as. 


£+io.^ 
^0  > 


(1) 


where 

M  is  the  mean  wind  (in  imits  m/ s). 


u*=jv 


du  . 
dz 


is  the  friction  velocity  (related  to  surface  shearing  stress). 


V  is  kinematic  viscosity  (in  units  m 
z  is  height  above  ground  level. 
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and 

Zq  is  the  roughness  (length)  of  the  surface  boundary. 

Surface  (aerodynamic)  roughness  of  the  landscape  affects  airflow  close  to 
the  ground.  The  surface  roughness  length  is  a  quantity  that  varies  in  the 
following  manner:  Zq  is  m  over  snow,  sand,  dry  lakebeds,  or  concrete, 
~10“^m  over  soils,  and  ~0.1  m  over  farmlands,  tall  grass,  or  shrubs  (e.g., 
Oke,  1978;  Pielke,  1984).  Typical  values  for  these  data  also  can  be  found  in  a 
table  compiled  by  Rachele  and  Tunick  (1994).  Knowing  z  and  Zq,  and  hav¬ 
ing  determined  from  the  sensor  retrieved  data  (as  explained  above),  we  can 
estimate  a  value  for  u*,  the  friction  velocity.  Then  the  application  of  Eq.  (1) 
allows  us  to  extrapolate  wind  velocity  profiles  (~15  m)  upward  through  the 
sensor  (fig.  2).  In  addition,  we  have  calculated  a  useful  turbulence  param¬ 
eter,  the  square  of  which  (in  simplified  form)  can  be  equated  to  the  surface 

momentum  flux,  i.e.,  u}  =  -u'  w'  is  the  rate  that  momentum  is  transferred 
downward  due  to  surface  friction.  Likewise,  the  square  of  the  friction  ve¬ 
locity  provides  an  estimate  of  the  surface  shearing  stress  (per  unit  mass  per 
unit  volume  of  air),  i.e.,  the  drag  force.  This  parameter  is  an  important  ele¬ 
ment  in  determining  the  surface  layer  budget  of  turbulent  kinetic  energy. 

Alternately  we  can  estimate  surface  layer  wind  velocity  profiles  using  a 
power-law  exponent  model,  such  as  described  by  Touma  (1977)  and  Irwin 
(1979),  i.e. 


m(z)  = 


f  y 


(2) 


where  is  the  mean  wind  velocity  at  reference  height,  and  p  is  specified 
as  a  function  of  static  (thermal)  stability  and  surface  roughness.  These  au¬ 
thors  have  provided  tables  to  give  examples  of  predicted  and  observed  val¬ 
ues  for  the  power-law  exponent,  p,  for  various  sites. 


Figure  2.  Schematic 
showing  wind  velocity 
profile  information 
provided  through 
application  of  a  simple 
log-law  model. 
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For  iiow,  we  present  the  simple  log-law  and  power-law  formulations  de¬ 
scribed  above  to  extrapolate  wind  profile  information  upward  through  the 
sensor.  In  future  algorithms  we  will  improve  upon  these  using  more  so¬ 
phisticated  wind  profile  models  for  the  lower-boundary  layer  for  day  and 
nighttime  conditions,  such  as  those  described  by  Santoso  and  Stull  (1998, 
2001). 

To  complete  the  algorithm,  we  the  calculate  the  surface  layer  vertical  tem¬ 
perature  profile  as  described  by  Hess  (1979),  i.e., 

f{z)  =  %  +  T{z-z,)  ,  (3) 

where  f(z)  is  the  mean  temperature  profile,  %  is  temperature  at  the  refer¬ 
ence  height,  and  T  is  the  adiabatic  lapse  rate,  i.e.,  F  =  -0.0098  °C/m.  In  the 
absence  of  buoyancy,  i.e.,  neutral  (windy)  conditions,  temperatures  in  the 
surface  layer  often  vary  linearly  with  height. 


4.  Summary  and  Conclusions 


A  precision  thermal  anemometer  and  meteorological  sensor  is  being  evalu¬ 
ated  for  implementation  onto  one  of  ARL's  all-terrain,  prototype,  robotics 
vehicles.  Compact  packing  is  important,  so  that  in  application,  this  meteo¬ 
rological  sensor  system  can  be  used  to  retrieve  and  construct  local  wind 
field,  temperature,  and  turbulence  information  alongside  infrared,  video, 
and  acoustic  arrays.  It  is  expected  that  instantaneous  meteorological  meas¬ 
urements  (with  the  use  of  atmospheric  computer  models)  will  contribute 
useful  information  on  the  performance  of  future  battle  field  sensors  and 
communications  and  thus  affect  elements  of  the  Army's  command  and  con¬ 
trol  architectures. 

In  this  report,  therefore  we  presented  a  new  algorithm  to  calculate  microme- 
teorological  profile  and  turbulence  information,  which  includes  a  simple 
data  analysis  algorithm  for  calculating  the  small  interval  (<  4.5  min)  aver¬ 
ages  and  standard  deviations  from  measured  data.  In  addition,  we  described 
an  algorithm  that  extrapolates  wind  field  and  temperature  profile  informa¬ 
tion  upward  through  the  sensors  platform  and  provided  estimates  of  an 
important  turbulence  parameter  related  to  surface  wind  shear.  As  a  result, 
we  have  intiated  a  coupling  of  meterology  to  robotics  for  future  military 
sensor  and  systems  applications. 
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